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TEE BESSELIAN FUNCTION, 



BY PROF. ASAPH HALL. 



The investigation in the theory of elliptic motion of the relations be- 
tween the different anomalies, the expression of the radius-vector by 
means of series, and the expansion into series of various functions of the 
radius-vector and the anomalies, have led to interesting and valuable re- 
sults. The theorem of Lagrange has hero a direct application, but the 
method of determining the coefficients by means of definite integrals is 
more elegant. This method appears to have been first indicated by 
Poisson in the first edition of his Mecanique, but it was first completely 
worked out by Bessel in a memoir published among the memoirs of the 
Berlin Academy for 1816-17. 

Let s and g be the excentric and the mean anomalies, and e the excen 

tricity of the ellipse, then 

g = e — e sin e. 

Assume e — g = A x sin g + A% sin 2 g + -j- A { sin * g+ .... (1) 

If i and *' are two whole numbers, we shall have when they are unequal 

/cos i g cos V g. d g = 0, I sin * g sin *' g . d g = 0, 
and when * = i', 

cos 8 ig . d g=\ it, J sin 2 i . dg = \n. 
Multiplying equation (1) by sin i g . d g and integrating, we shall have 
Ai == — j (e — g) sin i g . d g. 

It is more convenient to use the variable c, and integrating by parts' 
since e — g is zero at the limits o and it, we have 

or since I cos i g . d g = 0, 

2 C* 2 /** 

A t = -i- I cos i g . d e = — I cos z (e — e sin e) . d e . . . . . (2) 

For the radius-vector we have a series proceeding according to the 
cosines of multiples of the mean anomaly, and the coefficients in this 
series can be derived from the value of A, by a simple differentiation 
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with respect to the excentricity. In the series for the equation of center 
the coefficients can also be derived from the value of A t , but by a more 
complicated relation. 

In order to show another occurrence of this definite integral let it be 
required to change a series proceeding by sines and cosines of multiples 
of c into one proceeding according to sines and cosines of multiples 
of g, a problem that occurs in the theory of perturbations. Here he have 
to solve the equations 

sin m e = 2 AT sin i g; cos m e = 2 Bf cos i g: 

As before we have 

Af — — I sin m e sin i g . d g; Bf = — I cos m e cos * g . d g. 

Integrating by parts and omitting terms that are zero at the limits we 
have 
Af = -. — I cos iq co&m e . a e; Bf = -. — I sin i a sin to e . d e. 

Putting for g its value, e — e sin e, and then changing the products of 
the cosines and sines to sums and differences we find 

Af = ^- ( cos [(*' -f- to) e — i e sin e] . d e 

-4- -1- I cos [ (i — m) e — i e sin si . d e, 

Bf = — I cos [(i — to) e — % e sin s] . d s 

1/TtJ a 

— -j— I cos [ (* + m) e — i e sin e] . d e. 

These integrals are of the same form as that in equation (2). Hence if 

we put 

1 C n 
J£= — I cos (i e — k sin e) . d e (3) 

TtJo 

we have 

Af = ™ (Jif m + JiT™); b?=™ wr - Jit**)- 

If therefore we have a table of the J function we can compute the co- 
efficients easily. This function occurs also in the solution of the partial 
difierential equations which are found in the theories of wave motion and 
of heat; and as it was first investigated and tabulated by Bessel it is 
called by the German mathematicians the Besselian function. More com- 
plete tables have been computed by Hansen, who has investigated this 
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function in his peculiar way, and who has brought out many curious prop- 
erties. 

We can find an equation between three successive values of the func- 
tion as follows: Let 

u = sin (*' e — k sin e), 
then 

k 
d u = i cos (* e — k sin e) d e — 5- cos [(*' +1) e — k sin e] d e 

— -5- cos [(« — 1) e — k sin e] . d e. 

Since u is zero at the limits and » we have by integrating this value 
of du, 

kj£-i-2iJ£ + kJZ+i=Q (4) 

This equation gives the value of the function from the values of the two 
lower orders, but it is not well adapted to numerical calculation, since it 
gives the value of a small quantity from the difference of two greater 
ones. It can however be easily transformed into a continued fraction well 
suited to such calculation. 
From equation (3) we have 

J£ = — I cos i e cos (k sin s) . d e -| I sin * e sin (k sin e) . d t. 

Tt J o 71 J o 

And hence 

1 C* 1 f* 

j£ = — I cos (k sin s) . d e = — 8 d e 

Tt J o Tt J o 

Y \ 1 /<: 2 sin 2 e . &* sin 4 e ¥ sin 6 e . \ 

\ ~T7¥~ + 1.2. 3.4 1.2.3.4.5.6 + " ' " j ? 

1 C n If" 

J£ = — I sin e sin (k sin e) . d e — — I d s 

Tt J o * Tt J o 

X \ k sin 2 e _ Fsin4£ 4. *" sin6£ _ I 

XjAsin e- T1I + 1.2.3.4.5 •"• T 

Integrating the terms of the series between the limits and Tt we have 

T-o _ -, ** , ** » , & 

J h — 1 — -z -+■ — + 



2 2 ' (2.4) 2 (2.4.6) 1 ' (2.4.6.8)* 

jri _ 3* __ 4#> , 6& 5 __ 8ff , 

4 ~ 2 2 (2.4) ! + (2.4.6) 3 (2.4.6.8) 2 ~ i ~"" 

These values and equation (4) enable us to compute all values of the 
function. An elegant derivation of the general series for Jl will be found 
in Schlomilch's Compendium of the Higher Analysis, Vol. 2, p. 157. 
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This derivation depends on the following formula given by Jacobi, 
(Orelle, Vol. 15): 

fV(cos t) cos n t . d t = ,-s-g — ^ jr . Cf* (cos t) sin 2n t.dt, 

where (n) denotes repeated differentiations. 



DISCUSSION OF AN EXPONENTIAL CURVE. 




BY IKV1NG P. CHUEOH, B. C. E., NEWBURGH, N. Y. 

Construction by Points. — Let C = o be the base of any logarithmic 
system. On Cand G V, perpendicular to C, lay off O D = u, and 
C E = v, such that u = log v in the system whose base 
| is c. The intersection of B and D P, (parallel to 
O) fixes P, a point of the curve. 

From the construction we can immediately derive the 
I polar equation of the curve : 

q p = s _ _JfL = lo g^ = log • [g tan . ff| 
sin ^ siu ^ sin ^ 

_ log . o -\- log tan <p „, 

sin ^ 

Substituting from y = s sin </>, and tan ^ = -^-, we have 

OS 

for rectangular co-ordinates, 

y = log . d + log (JJ.) = log . | ^ | or <f> = ^, 

whence we have the final equation of the curve, 

<* = ^; (2) 

and also successively, 

p y = o^(l-ylo), 

(J T 

and -=—5 = c x ~ v I o (y I c — 2); here I e = Nap. log . c. 

By reference to the equation (2) we see that if c > unity the curve has 
the form indicated by our first diagram, which will give 
a maximum for x; also if c = unity, (2) becomes x = 
y, the equation of a straight line bisecting the first angle, 
while if g < unity, the curve passes to the other side 
of the line just mentioned, and is tangent to Y at y — 
— oo, thus, and gives a minimum for x. 




